Conformational Properties of an Adsorbed Charged Polymer 
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The behavior of a strongly charged polymer adsorbed on an oppositely charged surface of low- 
dielectric constant is formulated by the functional integral method. By separating the translational, 
conformational, and fluctuational degrees of freedom, the scaling behaviors for both the height of 
the polymer and the thickness of the diffusion layer are determined. Unlike the results predicted by 
scaling theory, we identified the continuous crossover from the weak compression to the compression 
regime. All the analytical results are found to be consistent with Monte-Carlo simulations. Finally, 
an alternative (operational) definition of a charged polymer adsorption is proposed. 
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Charged polymer (polyelectrolyte) adsorption on 
charged surface remains an interesting and important 
problem due to its influence to material science [jj], col- 
loidal science and biological science Hard sub- 
strates, and soft surfactant layers at interfaces can also 
be charged, due to the dissociation of ionic groups on 
the surfaces. Because the electrostatic force is strong 
and long-ranged, the electrostatic interaction between a 
charged polymer and a charged surface usually dominates 
over other non-electrostatic ones. 

The problem of charged polymer adsorption on 
charged surface can be studied by many approaches Q. 
By replacing the counterion effect by the Debye-Hiickel 
potential within the linear mean-field theory, one solve 
the Edwards equation 0, 0, • One can also solve both 
the Edwards equation and the Poisson-Boltzmann equa- 
tion self-consistently ^ Jii] at nonlinear mean-field 
level in which the effective screening length near the 
charged surface may not be equal to the bulk one. Scaling 
theory were also applied to the problem Even more, 
the effect of attractive image forces from high-dielectric 
substrate 0, lisl , and repulsive image forces from 
low-dielectric substrate [l^ [l^ [13, llal were also investi- 
gated by analytical methods or Monte-Carlo (MC) sim- 
ulation. 

However, the analytical approaches involving the Ed- 
wards equation usually impose zero monomer density at 
the charged surface in which the electrostatic bound- 
ary condition cannot be faithfully respected. It is only 
for the case of charged polymer adsorption on the high- 
dielectric substrate studied by Cheng et al. that 
the surface monomer density is properly treated. The 
surface monomer density follows a linear relation with 
surface charge density at Debye-Hiickel level. It indicates 
that the charged polymer is fully compressed on the high- 
dielectric substrate without any conformational change. 
For the low-dielectric substrate, due to the repulsive im- 
age forces, the polymer is not necessarily compressed on 



the substrate. Instead, the conformational degree of free- 
dom plays an important role on the adsorption behavior. 

In this paper, we study the conformational properties 
of charged polymer adsorbed on the low-dielectric sub- 
strate at Debye-Hiickel level by both the functional in- 
tegral methods and MC simulation. It is found that the 
usual Edwards equation is no longer valid to describe 
the non-Gaussian feature of polymer conformation. A 
new formulation by the functional integral method is pro- 
posed and compared with simulation results. Finally we 
give an operational definition of charged polymer adsorp- 
tion. 

A charged polymer carrying positive charges is im- 
mersed in a medium (z > 0) of dielectric constant e. 
At z = there is an impenetrable surface. Below the 
surface (z < 0), there is the substrate of low dielectric 
constant e' < e. Just above the substrate, there is an 
uniform surface charge density ct < 0. The adsorbed 
charged polymer always stays above the surface charge 
layer. Denote the charge on a polymer segment ds by 
qds, the Hamiltonian is 

T l-N rN ( K|f=-(s)-r(s')l 

n = -l ds I ds'lT^^ — _^ + r'(2-5,v) 
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where s is the variable to parametrize the chain and 
the Debye screening length. r(s) = {x{s),y{s), z{s)), 
r'(s') = {x{s'),y{s'),—z{s'j) are the positions of the 
monomers and their electrostatic images, respectively, 
r = gVe, F' = F(e - e')/{e + e') > 0, and h = 
4:TTq\a\/{e' + e) > arc the coupling parameters gov- 
erning the strengths of Coulomb interactions among the 
monomers themselves, between the polymer and its im- 
age, and between the polymer and the charged surface, 
respectively. Note that the above Hamiltonian is not ex- 
act even at Debye-Hiickel level. In particular, the longi- 
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tudinal interaction decays algebraically rather than ex- 
ponentially (17.] . However, the conformational properties 
related to the adsorption behavior will not be affected. 
We shall focus on the case of a charged polymer adsorp- 
tion in a low ionic strength medium. 

The continuum Hamiltonian in Eq. Q is discretized to 
perform MC simulation. The continuous curve r(s) is 
replaced by a chain of beads fi {i = 1, . . . , N) with hard- 
core excluded volume of finite radius a. The length and 
energy units are 2a and q^/2ea, respectively. Runs up to 
10^ MC steps and up to TV = 120 are performed. 

The partition function of the system is 
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where r\\{s) — (x(s),j/(s)) is the xy-plane projection of 
the curve r(s). While the charged polymer is adsorbed, 
lA^il <^ |Ar||j|, and note that f\\{s) should describe a 2D 
polymer conformation. Hence we approximate 
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wherefj_(s) = (|r|| (s)|, z(s)) is the side-view of f(s) along 
the curve f\\ (s). Note that the coefficients of the entropy 
terms of r||(s) and f±{s) are — 1/a^ and — l/2a^, respec- 
tively, which are different from that of r{s), — 3/2a^. 

For the case of charged polymer adsorption, the self- 
electrostatic interaction takes almost no effect in r±{s) 
since \r{s) — r{s')\ ~ 1*^1(5) ~ The repulsion from 

the images of the monomers can be effectively approxi- 
mated by the interaction between each monomer and its 
image only. The residual repulsion is absorbed by renor- 
malizing T' . Then the partition function becomes 
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The system is decoupled into two independent degrees of 
freedom, f\\{s) and rj_(s). Since the above functional in- 
tegral with respect to r|| (s) does not affect the adsorption 
behavior, we investigate only the conformational proper- 
ties of r±{s) in the following. 
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FIG. 1: Schematic diagram for the conformation of an ad- 
sorbed charged polymer. The degrees of freedom (entropies) 
of the polymer consists of three parts, the translation {r±c), 
the conformation (t(s)), and the local fluctuation {Sf±{s)). 
The solid line represents the polymer orientation. The local 
fluctuation lies within the blobs, (a) The weakly compressed 
polymer (onset of adsorption) conformation in general, (b) 
The weakly compressed polymer conformation in our analyti- 
cal approximation, (c) The compressed polymer (adsorption) 
conformation in which the conformational degree of freedom 
t{s) vanishes. 

Because of the repulsive image force from the low- 
dielectric substrate, the charged polymer may be at weak 
compression or compression in which their schematic di- 
agrams are shown in Fig^ and respectively. The 
terminology of weak compression (onset of adsorption) 
and compression (adsorption) are borrowed from Borisov 
et al. 0, for grafted polymers. 

In order to distinguish between the weak compression 
and the compression in our formulation, and note that a 
slowly varying orientation of polymer conformation un- 
der weak compression, we decompose 



fj_{s) = f_i_c + t{s) + 6f^{s) 



(5) 



where rj^c — jj Jq dsf±{s) is the position of the center 
of mass, and t{s) is the orientation vector of the charged 
polymer. We also restrict 



r|| (s' 



t{s)-Sr±{s) = 
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so that 5'r±{s) represents the local fluctutation along t{s). 

The adsorbed polymer is now characterized by trans- 
lational {r±c)i conformational (i(s)), and local fluctuta- 
tional {5'r±{s)) degrees of freedom. Under the compres- 
sion regime, t{s) vanishes. 

In general, it is hard to compute the effect from t{s). 
For simplicity but still capturing the qualitative picture 
of the weak compression as shown in Fig^, we further 
make an approximation that 
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and repeat for a period of 21 /a. Its schematic diagram 
is shown in FiglTJj. Substituting Eqs.|(Sl)-(|7I) into Eq.lgl), 
and at low-salt limit, we get 
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where the integral of (s) is absorbed into the normaliza- 
tion constant TV. Expand the following integral around 
small Sr±{s) up to quadratic order. 
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and then integrate out the variable t(s) under the condi- 
tion that |t • z| <C 1, the partition function becomes 



Z = N- 



/ dzcZc exp[-NP{hzc + - — )] 
Jo 4zc 
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Note that I is related to t{s) via Eq.{7|), and will be inte- 
grated out inside the functional integral of t(s). 

Without the effect from t{s) and (5rj_(s) which ex- 
pressed in the effective potential of 6z{s), the poly- 
mer acts as a rigid rod. Its equilibrium height is at 
Zc = ^JV l^h. The ensemble average 
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The entropic force points upward (downward) when the 
height of center of mass of the polymer is lower (higher) 
than v/ryiX. If F' = (same dielectric constants), the 
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FIG. 2: (a) Simulation results for the equilibrium height of 
the polymer, z'^^ (in units of 2a), as a function of surface 
charge density, \cj\ (in units of q/4,a^), in logarithmic scale at 
e'/e = 0.01 and = 25. It shows the scaling exponents 
decreases continuously from -0.96 (weak compression) to 
(compression) with the surface charge density. The scaling 
exponent of -0.5 indicating the crossover (analytically pre- 
dicted) occurs at \a\ ~ 0.11. (b) The thickness of the diffu- 
sion layer, Ds^, as a function of surface charge density, \a\, 
in logarithmic scale. It shows the scaling exponents are -0.67 
and -1.00 for weak compression and compression, respectively. 
The crossover occurs at \a\ ~ 0.11. 



effective potential for 6z{s) is linear rather than the har- 
monic. The result for the case of low-dielectric substrate 
cannot be analytically continued to the case of same di- 
electric constants. If F' < (high dielectric substrate), 
the system is unstable. It implies that the decomposition 
in Eq.|Hl) is inadequate in high-dielectric case. 

Hence, after integrating out the fluctuation variable 
5z{s) under the ground state dominance (large- iV limit), 
there leaves only the variable Zc in the partition function 
which determine the effective probability density distri- 
bution for the height of the center of mass. 
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up to a normalization constant. The new equilibrium 
including the effect from conformational changes is cal- 
culated by "force balance", dz^ \ogp{zc) — 0, which gives 



Kzl'^f^-A 
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at large- iV limit. For high enough surface charge density 
that z'^'^ is low, Ea. H13|l gives z"'' ^ |(t|°. The polymer is 
compressed in which the center of mass is independent 
of the surface charge density. When the surface charge 
density is lowered such that z^"^ is high, Ea. lfT!^ reduces 
to the scaling z°'^ ~ jcrl^^/^. If the surface charge den- 
sity is further lowered so that the polymer basically be- 
haves as a colloid (undeformed state), Ea. (|12|) becomes 
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exp[— iV/3/iZc], and hence z®i ~ |cr| ^. It pre- 



dicts a continuous crossover from the compressed state 
to the weakly compressed state. It is different from that 
obtained by scahng analysis for grafted polymer, predict- 
ing a discontinuous jump Our analytical result is 
consistent with MC simulation as shown in Fig|2K- 




FIG. 3: Monte-Carlo results for the normalized monomer den- 
sity at the surface, pa, as a function of surface charge density, 
a (in units of q/Aa^), at e /e = 0.01 and k"^ = 25. The 
straight Une is linearly fit to the data of higher \a\. It shows 
the data starts to deviate from the linearity when o ~ —0.11. 
Inset: More results to cover a larger range of surface charge 
density a. It shows the data follows the linearity at higher 



Besides the position of the center of mass, we also cal- 
culate the thickness of the diffusion layer which is defined 
as the characteristic length scale of the exponential decay 
of monomer density. We first determine the saddle-point 
0* from Eg. lfTUI) (equivalent to integrating out the vari- 
able Zc at large- iV limit), which is given by 



h{zlf-^z:+T'{5z{s))=Q 



(14) 



The above equation is then solved self-consistently with 
Ea. (|ll|l in which Zc is replaced by z*. The solution is 

z: = yrV4^^|a|-i (15) 

The effective partition function for 5z[s) becomes 

' d6z{s)'^ ^ 



Z = j V[5z{s)]eMj^ '^^^'^ 
which gives the diffusion layer thickness 



ds 
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is independent of the polymer conformation, which is 
consistent with its definition expressed in Eqs.(|SJ-l(ni. 
However, simulation results in Figl^h show that the scal- 
ing exponent is -0.67, a quite large deviation from our 
analytical result, -0.375. The deviation may be due to 
the approximation of the effective potential up to the 
quadratic order only. As shown in Figl^t- and|2t>, both 
the simulation results of z"'^ and Dsz exhibit the crossover 
between the weak compression and compression regimes 
occuring at cr ~ —0.11. 

Finally we also examine the relation between the sur- 
face monomer density and surface charge density by MC 
simulation. FigEl shows the simulation data follows the 
linearity at high enough |(t|, and start to deviate from 
the linearity at a ~ —0.11. The linear relation implies 
the compression regime The deviation from linear- 
ity tells that the polymer starts to be weakly compressed, 
which is also consistent with the MC results of both z'^'^ 
and Dsz- Since it is hard to characterize the polymer 
conformation in MC simulation by the original definition 
as shown in Fig^ we would like to propose an alternative 
(operational) definition for charged polymer adsorption 
- the linearity between the surface monomer density and 
the surface charge density. 
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The scaling exponent does not depend on the surface 
charge density. It means that the local fluctuation 5z{s) 
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